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Abstract 

The scattering-matrix product SS^ of a weakly absorbing medium is related by a unitary transformation to the 
time-delay matrix without absorption. It follows from this relationship that the eigenvalues of SS^ for a weakly 
absorbing chaotic cavity are distributed according to a generalized Laguerre ensemble. 
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1. Problem 

The purpose of this note is to answer a question 
raised by Kogan, Mello, and Liqun [IjJ, concerning 
the statistics of the eigenvalues of the scattering- 
matrix product SS> for an absorbing optical cav- 
ity with chaotic dynamics. Without absorption 
the scattering matrix S is an N x N unitary ma- 
trix, hence SS^ is simply the unit matrix. With 
absorption the eigenvalues R%, R 2 , . . . Rn of SS T 
are real numbers between and 1. What is the 
probability distribution P({R n }) of these reflec- 
tion eigenvalues in an ensemble of chaotic cavi- 
ties? 

In principle, this problem can be solved by 
starting from the known distribution of S in the 
absence of absorption (which is Dyson's circu- 
lar ensemble 0), and incorporating the effects 
of absorption by a fictitious lead j|. What has 
been calculated in this way is the distribution 
P({R n }) for small N || and the density p(R) = 
(£n S ( R ~ R ™)) for lar S e N B- These results 



have a complicated form, quite unlike those famil- 
iar from the classical ensembles of random-matrix 
theory || . For example, in the presence of time- 
reversal symmetry the distribution for N = 2 is 
given by g 

p(r x ,r 2 ) = (i-i?ir 4 (i-i? 2 r 4 

x exp[- 7 (l - Ri)- 1 - 7 (1 - R2)- 1 ] 
x \Ri-R 2 \ [7 2 (1 -e 27 +7 + 7e 27 ) 



- 7 (.Ri +R 2 - 2)(| - |e 27 + 2 7 + 7 e 27 
(l-i? 1 )(l-i? 2 )(3-3c 27 + | 7 



| 7 e 27 + 3 7 2 + 7 3 ) 



where 7 is the ratio of the mean dwell time Td 
inside the cavity^] and the absorption time r a . 

The situation is simpler for reflection from an 
absorbing disordered waveguide. In the limit that 
the length of the waveguide goes to infinity, the 
distribution of the reflection eigenvalues becomes 
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2 The mean dwell time is related to the mean frequency 
interval A of the cavity modes by = 2n/NA, so that 
7 = 2w(T a NA)- 1 . The definition of 7 used in Rcf. || 
differs by a factor N. 
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that of the Laguerre ensemble, after a transforma- 
tion of variables from R n to \ n = i?„(l — i?„) _1 > 
0. The (unnormalized) distribution is given by 



^({A„}) oc J] \Xi - Xj 



i<3 



J]exp[-7(/37V + 2-/3)A fe ], 



(2) 



where now 7 = r s /r a contains the scattering time 
r s of the disorder. The integer j3 = 1(2) in 
the presence (absence) of time-reversal symme- 
try. The eigenvalue density is given by a sum over 
Laguerre polynomials, hence the name "Laguerre 
ensemble" §. 

Kogan, Mello, and Liqun used a maximum 
entropy assumption Q to argue that a chaotic 
cavity is also described by the Laguerre ensemble, 
but in the variables R n instead of the A„'s. Their 
maximum entropy distribution is 



P({Rn}) oc J] \Ri - Rjf Y[exp{-aR k ) 



(3) 



The coefficient a in the exponent is left undeter- 
mined. 3 Comparison with computer simulations 
gave good agreement for strong absorption, but 
not for weak absorption This is unfortunate 
since the weak- absorption regime (7 >C 1) is likely 
to be the most interesting for optical experiments. 
Although we know from the exact small- TV results 
H that no simple distribution exists in the entire 
range of 7, one might hope for a simple eigenvalue 
distribution for small 7. What is it? 

2. Solution 

Absorption with rate l/r a is equivalent to a 
shift in frequency to by an imaginary amount 
8lj = i/2r a . If we denote by S{u>) the scattering 
matrix with absorption and by So (u>) the scat- 
tering matrix without absorption, then S(u>) = 
Sq(uj + i/2r a ). For weak absorption we can ex- 
pand 



5 (w + i/2r a ) 



ATr, CLLU 



3 We have verified that the theory of Ref. jij agrees with 
Eq. (H) for strong absorption (7 2> 1), with coefficient 
a = fj/3N. 



So(w) 



1 - — Q(w) 
2r a 



(4) 



where Q — ~iSQdS /dui is the time-delay matrix 
||. Since 5o is unitary, Q is Hermitian. The 
eigenvalues of Q , the delay times t\ , t% , . . . tn , are 
real positive numbers. Eq. (^J) implies that, for 
weak absorption, 



5(w)5 f (w) = S„(w) 



(5) 



We conclude that the matrix product SS^ is 
related to the time-delay matrix Q by a unitary 
transformation. This relationship is a generaliza- 
tion to N > 1 of the result of Ramakrishna and 



Kumar |10| for N — 1 (when the unitary trans- 
formation becomes a simple identity). Because a 
unitary transformation leaves the eigenvalues un- 
changed, one has R n = 1 — r n /r a , or equivalently, 
A„ = r a /r„ (since A n — ► (1 - Rn)^ 1 for weak 
absorption) . 

The probability distribution of the delay times 
in a chaotic cavity has recently been calculated, 
first for N = 1 and later for any N Jl|. 

The corresponding distribution of the reflection 
eigenvalues for weak absorption is a generalized 
Laguerre ensemble in the variables A ra , 



P({A n })ocn^-A J -|' 3 

i<j 

xnAr /2 exp[-i 7 /^A fc ] 



(G) 



The eigenvalue density is given in terms of gen- 
eralized Laguerre polynomials, hence the name. 
The corresponding distribution of the reflection 
eigenvalues is 



P({Rn}) oc J] \Ri - Rj 



\0 



exp[-± 7j 3iV(l - Rk)- 1 } 

(1 - i? fc )2-/3+3/3JV/2 



. (1 



(7) 



The first moment of this distribution is 
N~ 1 (J2 n Rn) = 1—7, independent of 0. One 
can check that Eq. M) is the small-7 asymptote 
of the exact result M for N = 2, fi = 1. 
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In the case N = 1 of a single scattering channel 
the distribution (Q) reduces to 

(7/?/2) 1+ ' 3/2 



P(R) 



■(1-R)- 2 -^ 2 



(8) 



r(i + /3/2) 

x exp[-i 7 /3(l - R) 

including the normalization constant. Wc have 
plotted this function in Fig. 1 for 7 = 0.1 and (3 = 
1,2. It is totally different from the exponential 
distribution P(R) oc exp(-ai?) of Rcf. @. For 
comparison, we have also included in Fig. 1 the 
exact N — 1 result |l4| (which is known only for 
(3 = 2): 

P{R) = (1 - R)- 3 cxp[- 7 (l - R)- 1 } 
x [ 7 (e 7 -l) + (l+7-c^)(l-i?)]. (9) 

It is indeed close to the small-7 asymptote (ph. 




3. Conclusion 

Summarizing, the distribution of the reflection 
eigenvalues of a weakly absorbing chaotic cavity 
is the generalized Laguerre ensemble (^) in the 
parametrization A„ = R n (l — Rn) . The La- 
guerre ensemble (||) in the variables R n , follow- 
ing from the maximum entropy assumption [Q, 
is only valid for strong absorption. For inter- 
mediate absorption strengths the distribution is 
not of the form of the Laguerre ensemble in any 
parametrization, cf. Eq. (|l|). In contrast, the dis- 
tribution of a long disordered waveguide is the La- 
guerre ensemble (|^) for all absorption strengths. 

The relationship between the reflection eigen- 
values for weak absorption and the delay times 
implies that the delay times r„ for reflection from 
a disordered waveguide of infinite length are dis- 
tributed according to Eq. (^) if one substitutes 
7 A n — > t s /t„. The implications of this Laguerre 
ensemble for the delay times will be discussed 
elsewhere. 
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Figure 1. Probability distribution of the re- 
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that is coupled to the outside via a single-mode 
waveguide (inset). The solid curves have been 
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1 — 7 = 0.9). The dashed curve is the exact (3 = 2 
result from Eq. (^) . 
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